{(f^) for a scalar field in 2D black holes: a new uniform approximation 



V. Frolov/'^)* S. V. Sushkov/'')t and A. Zclnikov^'^''^)* 
Theoretical Physics Institute, Department of Physics, University of Alberta, Edmonton, AB, Canada T6G 2J1 
'^'^^ Department of Mathematics, Kazan State Pedagogical University, Mezhlauk st. 1, Kazan 420021, Russia 
P.N. Lebedev Physics Institute, Leninsky pr. 53, Moscow 117924, Russia 

We study nonconformal quantum scalar fields and averages of their local observables (such as 
in a spacetime of a 2-dimensional black hole. In order to get an analytical 
approximation for these expressions the WKB approximation is often used. We demonstrate that 
I at the horizon WKB approximation is violated for a nonconformal field, that is when the field mass 

■ or/and the parameter of non-minimal coupling do not vanish. We propose a new "uniform approx- 

' imation" which solves this problem. We use this approximation to obtain an improved analytical 

CN) I approximation for {<^^)"'" in the 2-dimensional black hole geometry. We compare the obtained 

^ . results with numerical calculations. 

D 



Q 
o 



I. INTRODUCTION 



Calculation of local observables which characterize vacuum polarization in the gravitational field of a black hole is 
an old problem. Main motivation is connected with study of back reaction effects and constructing a self-consistent 
^ ■ model of an evaporating black hole. A lot of information was obtained concerning properties of fluctuations {(f'^Y^" and 
vacuum stress-energy tensor (T^j/)""" in Schwarzschild and Reissner-Nordstrom geometries. A superscript ron indicates 
^ I ■ that we are dealing with physical finite quantities after a renormalization has been done. Combination of numerical 
] calculations [1-4] and different analytical approximations [4-10] usually gives quite good results. Nevertheless, there 
I . still remains one problem. Analytical 'approximation' does not work in the vicinity of a black hole horizon if the 
fS| ' spacetime is not Ricci flat and the quantum field is not conformally invariant. It gives logarithmically divergent result 
, near the horizon in the Hartlc-Hawking state, where one expects finite and smooth behavior. In four dimensions the 
' general structure of this divergence reads (see e.g. [10]) 

{^T"-a,H\gu\), (r^.r ~^4^1n(|ff„|). (1) 

X 

Here ai and a2 are Schwinger-DeWitt coefficients. For massive scalar fields 

= (3) 



Analogously, in two dimensions it is 



5 J d^v ai 



(^2)"" ^ ln(|5.t|) , {T,.)- - H\9u\) ■ (2) 



ioU 



' 1 ,1 



where /i is the mass of the field and ^ is a parameter of non-minimal coupling with the curvature TZ. 

In order to obtain the analytical approximation Anderson, Hiscock, and Samuel [4] used WKB approximation for 
solutions of radial equations. Recently it was demonstrated that the breakdown of the analytical approximation (see 
discussion in Appendix A) is directly related to the breakdown of WKB approximation for so called zero modes, 
that is zero- frequency solutions [11,12]. By using improved approximation for zero modes it is possible to correct 
the analytical approximation near the horizon and to make all calculated quantities finite. The concrete methods 
proposed in [11,12] are based on the approximations to zero modes valid only in the vicinity of the horizon. 
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The aim of this paper is to develop a new uniform approximation scheme for calculations of zero-modes contribution 

to and (T^i/)"'". Although the method is applicable in arbitrary dimensions, here, for simplicity, we restrict 

ourselves by studying {(p'^)"'" near two-dimensional black holes. We consider different interesting examples of 2D black 
hole metrics, and keep the mass of a scalar field m and the parameter of non-minimal coupling ^ arbitrary. 



II. GREEN'S FUNCTION 



Consider a massive quantum scalar field with an arbitrary curvature coupling in a spacetime of the 2-dimensional 
black hole. The metric of the static 2-dimensional black hole can be written in the form 



F ' 



(5) 



where the function F(r) should possess the following properties: (i) F(r) vanishes at the horizon r = tq and (ii) F(r) 
tends to 1 at r — > oo. It is convenient to rewrite the metric (5) in the dimensionless form as follows 



f 



fix) = F{r) 



(6) 



where i = t/ro and x = [r — ro)/ro, so that the horizon is located at a; = 0. We shall also use the dimensionless mass 
m = iiro and curvature R = TZr^. By making Wick's rotation f — > zr in the metric (6) one gets the Euclidean metric 
of the form 



ds 



(7) 



It is well known that the property of periodicity of Green functions in the Euclidean time coordinate r with the period 
(3 = T~^, corresponds to quantum systems at finite temperature T. We consider the case of nonzero temperature T 
and nonzero surface gravity k = f!^/2\^_^. 

The Euclidean Green function GEiX,X') is a solution of the equation 



[□b - - ^i?] Ge{x,t;x',t') = -d{x,T;x' ,t') 
which in the metric (7) has the following form 



dr"^ dx y dx 



Ge{x, t; x', t') = -6{t - r') 6{x - x') . 



This equation allows a separation of variables, so that we can write 

oo 

Ge{X, X') = TQoix, x') + 2T ^ cos K(r - r')) Qnix, x'), 



n=l 



where a;„ = 27rnT and Qn are "radial" Green functions which are the solutions of the 1-dimensional problem 



Qn{x,x') = -6{x - x') . 



The Green functions Gn can be written in the form 

gn{x,x') = S^^\x<)s^^\x>) , 



(8) 



(9) 



(10) 



(11) 



(12) 



where = max(a;, x') and x^ = min(a;, x'). Here the radial modes Sn^ and Sn^ are two linear independent solutions 
of the homogeneous equation of motion 



s{x) = . 



(13) 



We choose si^^ to be regular at infinity and Sn^ to be regular at the horizon. The solutions Sn^ and Sn^ are normalized 
by the condition 
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The coincidence limit of the Euchdean Green function gives the unrenormalized expression for (<^^), which is 
ultraviolet divergent, 

oo 

unren 

= Ge{x, t; X, t') = TGoix, a;) + 2r ^ cos(a;„(r - r')) ^„(a;,x). (15) 

n=l 

(For convenience the points are split in time direction only, so x' = x.) Subtraction of ultraviolet divergencies gives 
the renormalized (</?^) 

(^'>re„ = ((^'>u„re„ " (V''>Ds) > (16) 

T — >T 

where {<^^)ds is the DeWitt-Schwinger counterterm which in two dimensions has the following form 

('^%s = -ian4^+2C). (17) 



Here a is equal to one half the square of the distance between the points x and x' along the shortest geodesic connecting 
them, C is the Euler's constant. For a massive fields the constant fi is equal to its mass. For a massless field it is an 
arbitrary mass scale parameter. A particular choice of the value of corresponds to a finite renormalization of the 
coefficients of terms in the gravitational Lagrangian. 

The stress-energy tensor can be obtained by applying a certain differential operator to the Green function: 

(T^,) = lim D^,Ge{X, X'). (18) 

In order to calculate {^p^) and (T^i/) one needs to solve the radial equation of motion (13) and to find the radial 
Green functions Qn{x,x). In a general case it is impossible to find exact solutions of this equation, and so one oftes 
uses various approximate methods. In case the scalar field has a nonzero mass,-'^ one of the standard approaches is the 
WKB approximation. For example, recently Anderson, Hiscock and Samuel have derived the analytic approximation 
for {(f^) and (T^i/) for a scalar field in the general static spherically symmetric spacctime using the WKB expansion 
for radial modes [4] . However, the serious obstacle for applying this method in black hole spacetimes is that the WKB 
approximation even in case m > breaks down at the horizon for low-frequency modes. 

In the section IV we discuss this problem in more detail. 

III. EXACTLY SOLVABLE MODELS 

Let us consider the radial equation of motion (13). Let us introduce new variables 

s{x) = ^-^,, (19) 



f{x) y ' yV2 

and rewrite this equation in the form 



u{v) = . (20) 



Here T{y) = f{x) and 



R--^-y^ + ^ ■ (21) 



The prime means the derivative with respect to y. 

First of all, let us discuss cases when the equation (20) can be solved exactly. 



^Generally speaking, the mass should be large enough. If the mass is small or vanishes the WKB approximation cannot be 
applied for modes with small numbers n, including the zero mode, n = 0. 
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A. Conformally invariant scalar field 



Consider the 2d conformally invariant scalar field, i.e., assume to = and ^ = 0. In this case the equation (20) 
takes the form 



u = 0. 



Two independent solutions of this equation are 



,,(1) -,,1/2 „,(2) _ _ 1/2 , 



}_ ^0,^+1/2 (2) _ 1 j.-a;„+l/2 

y ' "n>o — /?r— » 



/2a;, 



'2a;n 



(22) 

(23) 
(24) 



B. 2D dilatonic black hole 



Consider a scalar field in the 2D dilatonic black hole spacetime with (see [13]) 

fix) = 1 - e-^ 



Using the coordinate y we find 



and 



The equation (20) now reads 



R 



1 + 2/' 



1 + 2/ 



d^u f oj„ - 



,2 1 



4 



+ 



m 



+ 



u = 0. 



(Iy2 y y2 y{l + y) y{l + y)^. 
This equation can be solved in terms of hypergeometric functions [13]. Its two independent solutions are 



(25) 



(26) 



(27) 



(28) 



„(i) = i/2-n„ / _y_ 

" ^ Vi + y 

„(2) ^ i/2-n„ ^ 



F(a,6; 2a;„ + 1; y^), 

F(a,6;217„ + 1;-^), 
1 + 2/ 



The Wronskian of these solutions is 

l^[«W,42)]^^„ = 4i) 

where 



1) dul^^ du^n^ (1) _ r(2a;„ + l)r(20„ + 1) 



dy dy " 
0„ = -v/o^JTm^, 



r(a)r(6) 



(29) 



(30) 



and 
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IV. THE WKB APPROXIMATION FOR RADIAL MODES 



In the cases when the analytical solution is not known WKB approximation is often used. Let us write the radial 
modes Sn^ and Sn^ in the WKB form 

^vrnM-fm'-')- '''' 

where n{x) is a new unknown function. Note that these solutions are properly normalized to satisfy the condition 
(14). Substituting (31), (32) into (13) gives the equation for Q{x): 

02 2 ^ 2,^.pf^ l/'f^'' 1.W 3/2l7'2 

Q =a;„ + m/ + ^i?/+-^ + -^--^. (33) 
This equation can be solved iteratively with the zeroth-order solution chosen as 

f^W = 17o=K + mV]'^'. (34) 

The first-order solution is 

o(i) - o,o-o,?^/,l /'^o , 1 ff'^'o 3 

n^> _no + n,-no + — + --^ + --^---^ . (35) 

Analogously, the fc'th-order solution can be found 

n^''^ = no + cii + ... + nk, (36) 

The WKB approximation is based on an assumption that 

f^o > f^i > ••• > fifc- (37) 

This ensures the convergence of series X^^q and allows one to break the series in order to construct an approximate 
solution for O of the form (36). Substituting fl^''^ into (31), (32) gives the WKB approximation of fc'th order for the 
radial modes Sn''^^ ■ By and substituting Sn'^^^^^ into (39) and then into (15) one obtains the WKB approximation 

for (V3^)unren: 

OO 

(^'(a;))unren,WKB = TC™(^, ^) + 2T ^ cos(a;„(T - t')) (38) 



with 



arB(a:,x)=.WWKB(^),(2)WKB(^)_ (39) 



Let us discuss the validity of this approximation for zero-modes, i.e., when n = 0. For n = one has flo = mf^/"^, 
i.e., r^o ~ Then Vti ^ f'"^ f~^^'^ (sec Eq. (35)). Note that at the horizon the function / vanishes while /' and k 

remain finite. Thus in the vicinity of the horizon fig <C fii. Moreover, one can show that near the horizon inequalities 
fio ^ f^i <C . . . <C 0/j are fulfilled instead of (37). This means that the series Y^'j^=o^k diverges, and the WKB 
expansion breaks down. Thus, the WKB method can not be used for constructing Gn=o near the horizon. 

V. NEW UNIFORM APPROXIMATION FOR THE RADIAL MODES 

A. General discussion 

In this section we propose a new approximation which can be used instead of the WKB approximation. 
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We choose the time T normaUzation so that the function J^{y) tends to 1 at y — > oo. We assume also that the 
surface gravity does not vanish, so that J^{y) ~ j/ at y — > 0. Let us write J^{y) in the form 

^fe) = |M. (4„, 

Here c > is a positive constant, and h{y) is a positive function depending on c and having the asymptotics 

h{y)\y^0 = ho + O{y) , ho > 

h{y)\y^^ = l + 0{y-') . (41) 
Note that such the representation explicitly reflects the asymptotical properties of J^{y): 

J'{y)\y^0 = kiC-^y + O{y^), 

J^{y)\y^o, = l + 0{y-^). (42) 
The expression (21) for the scalar curvature now reads 

R = -ic + y)^-yic + y)^+yic + y)l^. (43) 



The equation (20) takes the form 



^V)u = 0, (44) 



dy^ V y{c + y) y{c + yY 

where 

In the special case c = 1 and h{y) = 1 we obtain the 2D dilatonic black hole with J^{y) = y/{l + y) and V = (see 
Eqs. (26), (28)). 

In a general case we will fix the constant c by the algebraic condition 

h'{0) = 0, (46) 

in order that to guarantee the regular behavior of V{y) near the horizon y = 0. This condition shows that the function 
h{y) near y = has the form: 



h{y)\y^0 = h + O{y^). (47) 



We rewrite the equation (44) as 



where we denote 



fu 
dy 



^-{Un+iV)u = Q, (48) 



Let us compare two terms Un and £y. In case ^ = the term disappears and the equation (48) takes the form 

(Pu 

-^2-UnU = 0. (50) 

Now assume that ^ does not vanish. Let us consider separately asymptotical regions near the horizon {y — > 0) and 

far from it [y oo). Using the asymptotical form (47) for h{y) wc can find that near the horizon V{y) = 0{y^) and 
Un{y) = C(y~^)) and hence in the limit y ^ the absolute value of the term £y is much less than the one of Un- Far 
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yu„(y) 




FIG. 1. The graphs of the function y^Un=o{y)', (a) ^ = 0; (b) ^ = 1/6; the graphs corresponds to the values m = 0, 0.1, 0.2, 
0.3, 0.4, 0.5, 0.6, 0.7, 0.8 from bottom to top. 




FIG. 2. (a) The graph of the function y'^V{y). (b) The ratio ^V/Un, n = 0, 1, 2, m = 0.1, ^ = 1/6; T = (47r)"^ 



from the horizon h{y) = 1 + 0{y~^), then V{y) — 0{y^^) and Un{y) = C>{y^^), and so we can again conclude that 
\CV\ < \Un\ in the limit y ^ oo. 

In order to estimate values of the terms and C/„ in the intermediate region < y < oo we consider the particular 
configuration: two dimensional "Schwarzschild" black hole with 



fix) ^ —-. (51) 
a; + 1 



By using the relation dx/ f{x) — dy/y we find 



y^xe"". (52) 
h as a function of x can be obtained by comparing two representations, (40) and (51), for /: 

yh 



Hv) 

so that 



c + y 



h(y{x)) =. ^±^. (53) 
(x + l)e^ ^ ' 

The constant c can be fixed by the condition (46), which gives c = 1/2. V and C/„ calculated for this metric are shown 
in the figures 1 and 2. In the figure 1 graphs for y^[/„=o are given for various values of m and ^. Note that in case 
m < 0.5 the function y^Un=o{y) is strictly negative, while for larger masses m > 0.5 it can change sign, so that there 
is a point > where ytUn=o{y*) = 0. Hereafter we will assume that to < 0.5. As for the n > modes one can 
easily see from Eq. (49) that the function y'^Un>o{y) is positive for arbitrary m > and ^ > 0. The graph of y^V{y) 
is given in the figure 2(a). Note that the function V{y) is completely determined by h(y) and does not depend on the 
parameters to and In the figure 2(b) the ratio £,V/Un is shown for n = 0, 1, 2, to = 0.1, and ^ = 1/6. One has 
\^V/Un\ < 1. For example, max(|^y/C/„|) « 0.032, 0.009, 0.0002 for n = 0, 1, 2 respectively, so that the higher the 
number of mode the less the ratio |^y/C/„|. 

Thus we may conclude that at least for 2D Schwarzschild metric the term £,V{y) is much less than Un{y) in the 
whole range of y. Using the smallness of in comparison with J7„ we can neglect this term in the equation (44). 
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B. The uniform approximation for the massless field modes 

Consider the massless case, m = 0, with arbitrary couphng ^. Neglecting the term in the equation (44) we 
obtain 



or, after rescaling y = cy, 



+ 



y{c + yf 



Lot, 



+ ■ 



u = 0, 



u = Q 



Comparing with Eqs. (28), (29) we find two independent solutions of the equation (54): 



,(1) 



,.(2) 



(c + y)-''^"F(a, 6; 2a;„ + 1; ^) 



j/i/2+-„(c + y)-2-nF(a,6;2a;„ + l; 



c + y 

c + y' 



with 



a = i+2w„ + iVl-4e, 6=^+2w„-iVl-4e 
These functions represent an approximate solution of the equation (44) in case m = 0. 



(54) 



(55) 



(56) 



C. The uniform approximation for the msissive field modes 

Consider the massive case. Neglecting the term in the equation (44) gives 



dy"^ 



'h{y) 



y{c + y) y{c + yy 



u = 0. 



Generally this equation cannot be solved exactly, and so let us construct an approximate solution. For this 
study asymptotical properties of the equation (57). First consider the region far from the horizon, y ^ oo. Nej 
terms of order y~^ and smaller in the equation (57) we get 



(57) 

aim we 
fleeting 



rf2 ^ cjl-\+m'^h ' 



u = 0. 



(58) 

Introducing in the equation (58) a new variable^ by da;* = dy/y and using the relation s(a;*) = y~^/'^u{y) we obtain 

(59) 



^ - nl{x.)s = 0, 



where 



An approximate solution of the equation (59) can be written in the well-known WKB form: 



^/2^; 



exp 



■\/2iln 



exp 



(60) 



(61) 



^The variable x* is known as a tortoise coordinate, dx/ f{x) = dy/y = dx». 
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Returning to the variables y and u„ we obtain the WKB solution of the equation (58): 



,.(1) 



exp 



J y . 



,.(2) 



exp 



■ / 

J y . 



(62) 



These functions give an approximate solution of the equation (57) in the region y ^ oo. 

Now consider the region near the horizon, y 0. Taking into account Eq.(47) and neglecting in the equation (57) 
terms of order y~^ and smaller we get 



dy' 



u = 0. 



(63) 



It is worth noting that the last equation does not contain the parameters m and ^. This means that any scalar field 
with arbitrary mass m and coupling ^ near a horizon behaves effectively like the conformal scalar field with m = 
and ^ = 0. Two independent solutions of the equation (63) are 



"n=0 ~ y ' 
.(1) 



^n>0 



(2) _ 

w„ + l/2 



-y 



■ y 



1/2 In 

^n>0 ' 



y 



a + 1/2 



/2uJn v'2w„ 
Respectively, for the radial modes Sn''^^ = y~^^'^Un''^^ expressed in the coordinate y [see Eq. (19)] we obtain 



*n=0 
*n>0 



5(2) 



Iny 

(2) 
'^rt>0 



Note that the modes sh, ' are normalized by the condition (14) which for the coordinate y reads 



" dy 



(2)a£«_ 
" dy 



1 

y' 



(64) 
(65) 



(66) 
(67) 



(68) 



Thus, we may summarize that a solution of the equation (57) should possess the asymptotical properties (62) and 
(66). We also assiimod that terms containing derivatives of h are small and can bo neglected. At last, it should be 
taken into account that in case h = const the equation (57) could be solved analytically with the solutions given by 
(29).^ Now combining these results we choose approximate solutions of the equation (57) as follows: 



1/2 



with 



,(2) 



y^ 

\/2n^ 



exp 



exp 



dy_ 
y 



y . 



(69) 



c + y, 



F(a,6; 2w„ + 1; ^— ) 
c + y 



■y 



F{a,b; 2f2„ + 1; ■ 



y 



(70) 



21 

It can be shown that the given functions uii ' possess all necessary properties. First, substituting the solutions (69) 
into the equation (57) one may see that they obey the equation up to terms containing derivatives of h. Second, using 
properties of the hypergeometric functions (see Ref. [15]) one can verify that the solutions (69) have the asymptotical 



^Morc exactly speaking, in order to obtain the solution of the equation (57) in the form (29) in case h = const one has to 
make rescaling = m^h and y = c~^y. 
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form (66) and (62) at j/ — > and y — > oo, respectively. Moreover, note that the functions (69) reproduce the exact 
sohition (29) in case h = const or/and m = 0. 

Finally, using the relation s(x) = y~^^'^u{y) and the formulas (69) gives an approximation for the radial modes 
s„ ' \ It is worth noting that this approximation has been derived for arbitrary mass m of the scalar field (including 
TO = 0) and arbitrary coupling ^. It is also important that the approximation works correctly both far from and near 
the horizon. For this reason we will call it as a uniform approximation. 



VI. EVALUATING OF {ip'^} IN THE UNIFORM APPROXIMATION 

Using the uniform approximation, we may obtain the radial Green functions Qn = Sn^ Sn^ and construct the 
approximate expression for (</'^)unren- 

oo 

(¥'')u„re„ = TG^^i^, x) + 2Tj2 cos(27rnTe) ^^C^- (71) 

n=l 

where the superscript 'ua' is used to denote the uniform approximation. 

We may simplify the resulting approximate expression for (</5^)unren account that the n > modes 

are satisfactory described by the WKB approximation. In this case we will use the uniform approximation in order 
to compute Gn=o only, and use the WKB approximation in order to compute the other radial Green functions Qn>o, 
so that 

oo 

(^')unren = TG^^, x)+2TJ2 cos{2nnTe) G^'^'^ix, x). (72) 

n=l 

The corresponding calculations for (T^) would require to use the uniform approximation forn = 1 as well. Subtracting 
(v^)ds from (v5^)un].en ^"^^ taking the limit e ^ gives the renormalized expression for (</5^): 

(^'>ren = (V''>0 + (^')l. (73) 

where we denote 

{v^)^=TG^\x,x), (74) 

and 



{if )i = lim 



2r^ cos{2T,nTe)G^'^''{x,x) - {^^) 



DS 



(75) 



Using the uniform approximation (69) for the radial modes s = y ^/^w and taking into account that 

ujo = 0, fio = mVh, Wo = — , 

r(a)r(6) 



we can write (<^^)o = TGq'^ = Tsf^^'s'^^ as follows: 

TT{a)T{b) ( 1 y-V^ 



r(2TOvW + i) vc + y 



Ffa,6;l;^')Ffa,6;2TOv^+l;^V (76) 
V c + yj V c + yj 



For a given spacetime the function h and the parameter c are known. Hence the expression (76) represents an 
analytical formula describing an approximation for {^p'^)^. 
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To evaluate {'p'^)i one may use the well-elaborated summation procedure dealing with the WKB approximation 
(see, e.g., [3,4], and also [14]). Combining Eqs.(33,39,75,A4), we obtain 



In ( / •^.^ 1 + 2C 
V2(27rr)2 



(77) 



where 



(27rT)2 



(78) 



and 



\ [(27rTn)2 + m^] 



2 fifc+V 



172' ^-1,2, 3,. 



(79) 



All the sums Sj in this formula converge everywhere and tends to zero at the horizon (/ 0). The logarithmic 
term combined with {^p'^)^, which also logarithmically diverges at the horizon, leads to (T — l/47r)ln(/) term (see 
Appendix). In the Hartlc-Hawking state T = 1/Att and, hence, total (i^^)r(,n finite, as it should be. 

Thus the uniform approximation solves the problem of finiteness of ((/?^) at the horizon by accurately taking into 
account of the zero mode contribution. Nevertheless finiteness itself on the horizon does not guarantee the correct 
value. So, we should compare our approximation with numerical calculations to get an idea how precise it is. The 
most important contribution comes from zero modes, this is why we present here plots for comparison of only their 
contribution. For higher modes the WKB approximation works very well (see, e.g., [4,14] One can see that the larger 



< 








0.95:^ 


1 


0.9:' 




0.85:^ 
0.8 + 
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FIG. 3. The graphs demonstrate the ratio of exact numerical calculations and the uniform approximation for {(p'^)q- The 
figures (a) and (b) correspond to ^ = and ^ — 1/6, and the curves a, b, c, d correspond to m = 0.05, 0.1, 0.2, 0.3, respectively. 



mass the more accurate becomes the uniform approximation. For m > 1 and ^ > 1/6 the uniform approximation 
gives practically exact value on the horizon and at infinity and the maximum deviation from exact function is less 
than 0.25%. For other parameters it is still very precise. This you can see in Fig. (3). 

Let us summarize the obtained results. We demonstrated that the WKB approximation traditionally used for 
obtaining different analytical approximations for local observables breaks down for low frequency modes. In the 
general case this results in the logarithmic divergences of these observables at the black hole horizon. The adopted 
approximations can be improved if one uses a more accurate expression for the low frequency modes. We propose a 
method, which we call a uniform approximation, which not only automatically gives a finite value for local observables, 
but also provides a good approximation of the observables uniformly in the interval from the horizon to infinity. Our 
concrete calculations were done for {(p'^) in 2D black hole spacetime, but the proposed method can be used for other 
local observables and for higher dimensional black holes. 
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APPENDIX A: REMARKS ON THE log-DIVERGENCE OF (93^) AND (T^^) NEAR THE HORIZON 

In this appendix we discuss the problem which accompanies many of approximate methods [see, for example [6,4,9]] 
derived for calculating expectation values, such as (<^^) and (T^v), on the black hole background. The essence of the 
problem is that the approximate expressions for these values turn out to be proportional to In / and hence diverge 
logarithmically near the horizon where / — > 0. It is important to stress that such the log-divergence takes place for 
both vacuum and thermal expectation values with any temperature, including the Hartle-Hawking state with the 
Hawking temperature. We will call this problem as the problem of log-divergence. 

As we will see the problem of log-divergence is tightly connected with accurate consideration and taking into account 
a contribution of zero modes into expectation values near the horizon. Really, as has been shown above the radial 
modes of a scalar field with arbitrary mass and coupling near the horizon have the asymptotical form (47). Then, the 

radial Green functions Qn = Sn"^ Sn"^ near the horizon are 

a„=o = -lnt/, (Al) 

Gn>o = -X — = -; — ^, (A2) 
and the unrenormalized expression for (<^^) reads 

, 9> 1 v"^ cos(27rnTe) _ 

iv Unren = -Tiny + — ^ ^, (A3) 



2tt 

n=l 



where e = r — r'. Using the formula 

cos(nKe) 1 



In 



1 

n=l 



= -ilnK2e' + 0(e'), (A4) 



_2(l-cos(Ke)) 

we obtain the asymptotical form for (<^^)unren near the horizon: 

(V5')u„ren = "Tiny - In K^e^ + Qie^), (A5) 
47r 

where k = 2ttT. To calculate a renormalized expression for {(p"^) one should subtract from ((p^)unrcn the DeWitt- 
Schwinger counerterms (17) and go to the limit e ^ 0: (</5^)].en = l™e^o (('^^)unren ~ (v^)ds)- Note that (</5^)ds 
— (47r)~^ Incr, and a /e^ in the limit e — > 0. Hence 

(¥'')ren ~ -^Iny + {in)-' In/ + 0{y'>). 

Taking into account that / « y at y ^ we rediscover the known result that (</'^)ren regular near the horizon 
provided the temperature has the Hawking value, T = (47r)"^, otherwise (v'^)ren diverges as Iny. 

It is worth noting that the term (Air)^' In / coming from the DeWitt-Schwinger counterterms (95^)153 is compensated 
by the term Tiny coming from the expression for the n = radial Green function Qo- Thus we see that the zero 
modes play an important role in the problem of log-divergence. 

To illustrate this statement we consider the approximation derived by Anderson, Hiscock and Samuel [4]. They 
obtained the approximate expressions for (9?^) and (T^^) [sec Eqs. (4.2), (4.3), (4.4)] which are proportional to In/ 
and hence diverge logarithmically near the horizon where / 0. The cause of this divergence can be easily explained 
now. Examining the procedure derived by Anderson, Hiscock and Samuel for computing (<^^) and (I),;/) for the 
nonzero temperature case one faces with a necessity to calculate the mode sums over n. Since the authors use for 
this aim an expansion in inverse powers of n, they encounter the problem of including into the calculation scheme 
a contribution of the n = modes. To resolve this problem the authors impose a lower limit cutoff in the mode 
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sums, which for the nonzero temperature case means merely discarding the n = contribution from the summation 

procedure (see a discussion in the section IV and the appendix E). But as was shown above, taking into account 
the n = modes is necessary for a compensation of the log-contribution coming from the counterterms (</'^)ds ^^'^ 
{Tij,v)bs', otherwise, one gets the problem of log-divergence. 
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